function meanM2 = model_CalMeanM2(model)
meanM2 = sum(model.meanAlpha2.*model.M2ModeWeighting(:,1)+model.meanBeta2.*model.M2ModeWeighting(:,2))/(2*model.SiteNum);